Interference and interaction effects in multilevel quantum dots 
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Using renormalization group techniques, we study spectral and transport properties of a spinless 
interacting quantum dot consisting of two levels coupled to metallic reservoirs. For strong Coulomb 
repulsion U and an applied Aharonov-Bohm phase <f>, we find a large direct tunnel splitting |A| ~ 
(r/-7r) | cos(0/2) | ln(U /cu c ) between the levels of the order of the level broadening V. As a consequence 
we discover a many-body resonance in the spectral density that can be measured via the absorption 
power. Furthermore, for (f> — tv, we show that the system can be tuned into an effective Anderson 
model with spin-dependent tunneling. 
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Introduction. Electronic transport through ultra-small 
quantum dots (QD), where the charging energy is the 
largest energy scale, has been studied extensively over 
the last few yearsE Due to the quantization of charge the 
transport is dominated by Coulomb blockade (CB). More 
recently experiments revealed that the transport can-be 
even more intriguing by measuring the Kondo effectcm as 
suggested in Ref.'s.cl 

The Kondo effect occurs for a dot with one low-lying 
spin-degenerate level. In this paper, we will study a dot 
consisting of two levels without spin or, equivalently, two 
dots in an Aharonov-Bohm (AB) geometry with one level 
per dot in the presence of an interdot Coulomb repulsion 
U . Such a system is of fundamental interest since the 
two possible paths through the dot (via level 1 or 2) can 
interfere with each other. The interference can be con- 
trolled by an AB flux and has attracted much interest 
due to the possibility of realizing AB interferometersu 
or use the coherent properties in connection with quan- 
tunKjeomputinga (for recent experimental realizations see 
Ref.LJ). Furthermore, in many recent experiments per- 
formed in the strong tunneling regime, the level broad- 
ening is large and transport is inevitably controlled by 
multilevel physics. 

The model. We consider a quantum dot consisting 
of two levels, labeled by j = 1,2. Via tunnel bar- 
riers the dot is connected to two electronic reservoirs 
r = L, R. The orbital index j is not conserved during 
tunneling and hence does not exist in the leads. The 
Hamiltonian is written as H = H^ot + H ICS + Ht, with 
Un]_n 2 , H les = Y,kr £ kral r akr, and 
+ H.C.). The tunnel matrix elements 



dot. 



are assumed to be real except for an AB phase, i.e., we 
attach a phase factor e 1 ^ to t\. The energy scale of the 



level broadening is defined by — 2n\t^\ po, where po is 
the density of states in the leads, which we assume to be 
independent of energy for the energy range of interest. 

We neglect spin (assuming a large Zeeman splitting) 
since the aim is at analyzing explicitly the physical effects 



arising from the tunneling-induced interference between 
the two levels. Since both levels overlap with the reser- 
voir states, there is an effective overlap matrix element 
—A/2, which, surprisingly, is shown to be zero for a non- 
interacting quantum dot, but for strong on-site Coulomb 
repulsion U >> \e\,V, scales like 



A 



v/ffrf 



L p i<f> 



In (U/u c ) . 



(1) 



Here, <j> is the AB phase, and uj c denotes a low-energy 
cutoff set by the maximum of the mean level position 
e = (ei+£2)/2, the mean level broadening r = (ri+r2)/2 
(with Tj = Tj 1 + rj"), the temperature T, or the bias 
voltage eV. The level splitting is given by 



8e= y^SiF 



|A| 2 , 



(2) 



where 5e = €2 — £i denotes the level spacing. Conse- 
quently, the tunnel splitting gives rise to an interference- 
and interaction-induced level repulsion, i.e., an effect not 
being considered in models with lewels labeled by a con- 
served quantum number (e.g. spinp or in the absence of 
interactions.u The energy scale of A is given by T and 
will influence the spectral properties as well as the con- 
ductance for low enough temperatures T < T. We em- 
phasize that this energy scale is well separated from the 
Kondo temperature Tk ~ VTU exp(-7re/r) (e -C — r), 
which is exponentially small and determines the crossover 
to the, occurrence of the Kondo effect for spin-degenerate 
levels! Most importantly, we will show in this paper that 
for low lying levels e — T (where the ground state is the 
singly occupied state), the effective level splitting shows 
up in a many-body resonance in the spectral density at 
the energy which e.g. can be measured by an absorp- 
tion experiment but influences also the temperature and 
flux dependence of the linear conductance. For <f> = -k 
and Tj* = Tj 1 , the tunnel splitting is zero, and the sys- 
tem is shown to be equivalent to an Anderson model with 
Zeeman splitting 5e. Thus, Kondo physics can be real- 
ized in a quantum dot without spin even if the quantum 
number labeling the levels is not conserved. 
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We note that multilevel ctata-in the presence of spin 
have been studied previouslyElHl However, Ref. ^studies 
the case of a conserved quantum number labeling the 
levels, and Refs. and |ll| consider the cases Se ^> T 
or tj> = ir, where the effect of the tunnel splitting A can 
be neglected. The same applies to the AB geometry of 
Rcf. |l2| where the interdot Coulomb repulsion is absent. 

Renormalization group study. An effective dot Hamil- 
tonian can easily be derived from perturbation theory 
or, equivalently, by integrating out the reservoir states 
by the renormalization group. The dot is characterized 
by four states |0), |1), |2), and |12), with energies E = 0, 
Ei = e\, E2 = £2, and E12 = e\ + e 2 + U. Intuitively, the 
hybridization with the reservoirs will lower the energies of 
all these states. For the singly occupied states, however, 
this is less pronounced because it costs a finite energy 
U to occupy the dot with a second electron. Therefore, 
the level positions ej — Ej — Eq will be renormalized 
upwards. Furthermore, a coupling between the levels is 
generated since tunneling events can shift the electrons 
between the two levels. For an electron starting in level 
1 there are two possibilities: either the electron first tun- 
nels out and hops into level 2 or an electron first hops 
into level 2 and then the electron tunnels out of level 1. 
In the latter case, the intermediate state is the doubly 
occupied state and, due to Fermi statistics, the matrix 
element gets an additional minus sign. Therefore, for 
reservoir electrons with an energy |efc| >> U, these two 
terms will cancel each other and, consequently, there is 
no direct coupling between the levels in the noninter- 
acting case. In contrast, for an interacting system, the 
doubly occupied state is suppressed, and there is a fi- 
nite coupling A between the two levels. We note that 
this mechanism does not work for levels characterized by 
spin since the two tunneling processes described above 
would also change the spin in the reservoirs and, there- 
fore, do not lead to a direct renormalization of the dot 
Hamiltonian. 

Using the real-time renormalization group (RG) of 
Ref. O for the forward propagator we find that energy 
scales lo c > U do not renormalize the states, i.e., we start 
the RG at ui c = min(D, U) where D is the bandwidth. In 
the basis of the three remaining states |0), |1), and |2), 
we obtain the flow equation (t c — 1/lu c ) 



dH, 



dot 



1 



dt r 



2n(t c - i0+) 



r! + r 2 
r x $ 
r 2 



(3) 



where $ = \/TfTf + ^T^T 



Neglecting level 



broadening, the solution of this equation gives an up- 
ward level shift E\ji — Eq = e-y^ + Ar 2 /i, with A = 
(1/271") \n(U/u) c ), and a coupling —A/2 = — $A leading 
to Eq. ([j]). As a consequence we get two effective levels 
at £1/2 = e + Ar =F Se/2, where the effective level splitting 
Se is given by Eq. While ei is quite close to the orig- 
inal level position, e 2 is strongly renormalized upwards. 
For ps T/2, Se < T and 4> < 1 the lower (upper) level 



is coupled strongly (weakly) to the reservoirs. For the 
following discussion we will usually assume equal cou- 
plings = T/2, i.e., \tj\ = tj = t and discuss the effect 
of asymmetries at the appropriate places. 

In the symmetric case we define y/2fi = c\ — (— l)*c 2 . 
For cj> = 0, only the /1 operator couples to the reser- 
voirs, whereas the level spacing Se controls the coupling 
between the /1 and / 2 level. The current operator in the 



H.c. 



right reservoir is given by Ir = ie\ 
We also note that for Se = and <f> — ir the conductance 
is exactly zero since the /1 (/ 2 ) level couples only to the 
left (right) reservoir. This is an effect of destructive inter- 
ference which interestingly persists also in the presence 
of interactions. 

Spectral density and absorption power. In Fig. [I] we 
show the spectral density of the /1 level for e > 0, where 
the ground state is given by the empty state. The results 
are obtained by using the full real-time renormalization 
group method of Ref. |l3| which is known to yield excellent 
results in the regime where charge fluctuations dominate. 
The two peaks in the spectral density correspond to the 
renormalized level positions and change qualitatively as 
function of temperature T and Se according to Eqs. (|l|) 
and (|^). The distance between the resonances saturates 
for Se < T at the energy scale A according to Eq. (||). 
In contrast, when e is below the Fermi level, the lower 
level is occupied and particle excitations lead to a broad 
shoulder in the spectral density at the effective spacing 
Se, see inset of Fig. |^ (an additional weak feature occurs 
at negative frequencies but this is masked by the broad 
resonance at ii). These results have been obtained by 
using Wilson's nopcerturbative numerical renormaliza- 
tion group (NRG)Hj which, up to some overbroadening 
effects at higher frequencies, gives very precise results for 
the spectral density near the Fermi level and for the posi- 
tions of all resonances. Since the location of the shoulder 
is not at the Fermi level, it is more suitable to test its 
position via the absorption power rather than the linear 
conductance. Therefore, we have shown in Fig. ^| the 
result for the spectral density of the transition operator 
c\c2 + c\c\. The peak position of the absorption power 
agrees very precisely with the position of the shoulder 
in the spectral density of the j\ level. We emphasize 
that the shoulder is absent without the tunnel splitting, 
i.e., it is a generic effect which will also be present in 
the asymmetric case Ti ^^2- In this case, however, the 
broadening of the shoulder (which is determined by T) 
will increase relative to its height (which is determined 
by A). 

Figure |^ shows the spectral density of the j\ level 
for different AB phases 4> (also obtained by NRG). For 
Se = 0, the position of the shoulder varies propor- 
tional to |A| - (r/7r)|cos(0/2)|ln([//|e|), according to 
Eq. (|l]) (for cf> = 0, the amplitude of the shoulder is 
zero since the /2 level is decoupled from the reservoirs). 
Furthermore, the resonances at finite frequency become 
more pronounced and, for cj> = tt, merge into a Kondo 
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resonance at the Fermi level with width given by the 
Kondo temperature Tk- This effect can easily be under- 
stood, since for (f> = tt, the tunneling Hamiltonian reads 
H T = tY,kj( b kj c j +H.c.) with \[2bu = a kR - (-l) l a kL . 
Hence, for this special case, the pseudo spin j is effec- 
tively a conserved quantum number and we obtain the 
Hamiltonian of the usual Anderson model, which^Jpr a 
low lying level e is equivalent to the Kondo modelEj 

We note that this realization of Kondo physics with- 
out explicit spin degrees of freedom is.-quite different 
from other realizations, where metallidl3 or two-level 
systemsEZl have been used. Furthermore, there are three 
experimentally tunable ways to destroy the Kondo res- 
onance. First, a finite level spacing 5e ^ acts like an 
effective Zeeman splitting. This splits the Kondo res- 
onance and decreases its height, see the left inset of 
Fig. [| Second, an AB phase away from <f) — it leads 
to an effective coupling A between the two levels. At 
A ~ Tk a phase transition will occur quite analog to the 
competition betweep-JlKKY and Kondo physics in two- 
impurity models Ej'Ej The same mechanism is induced 
by left /right asymmetries, i.e., for ^ T^. Third, for- 
given left/right symmetry but T\ ^ T2, we obtain an An- 
derson model with pseudo-spin-dependent tunneling ma- 
trix elements tj . As shown in the right inset of Fig. the 
Kondo resonance arising at $ = 7r is reduced and splits 
asymmetrically but is well defined even at T2/T1 w 2 (we 
note that the reduction is quite more pronounced for a 
finite level spacing). As a consequence, the Kondo reso- 
nance can be shifted away from the Fermi level by chang- 
ing the asymmetry of the tunneling matrix elements, an 
effect also seen in recent experiments J3 

Since the conductance is zero for 5e = and ef> = tt, 
the Kondo resonance will show up only weakly in the 
I(V) characteristics by changing the level spacing or the 
AB flux. However, the crossover to the Kondo effect can 
e.g. be measured by the absorption power. Alternatively, 
in an AB geometry with two dots and one level per dot 
we expect in equilibrium for Se — and 4> = tt a Kondo 
resonance in each dot separately. Their effect might be 
tested by measuring the conductance fluctuations of a 
parallel quantum wire lying very close to one dot. 

Linear conductance. Another fingerprint for the renor- 
malization of the energy levels due to the tunneling split- 
ting A is the measurement of the linear conductance. It is 
calculated by using the renormalized Hamiltonian on the 
forward and backward propagator according to Eq. (|^), 
including the level broadening (for ths-backward propa- 
gator we take the hermitian conjugated). This effective 
Hamiltonian is used as an input for the calculation of 
rates in lowest order in the tunneling coupling. 

Figure ^ shows the temperature dependence for Se = 
cf) = 0.113 At T = 0, the spectral density of the f\ level 
is a single Lorentzian with width Y centered at the level 
position e. The resonance at e+U is missing for e > —U /2 
since the fi level is decoupled from the system and is 
not occupied in the ground state; see also inset of Fig. 0. 



Thus, at zero temperature, the conductance is symmetric 
under a sign change of e, in contrast to the case for spin 
degenerate levels, where the Kondo effect enhances the 
conductance for negative e. For finite temperature, the 
ji level starts to become occupied and suppresses the 
conductance due to the Coulomb repulsion U. This effect 
is more pronounced for negative e and, therefore, the 
conductance shows a local maximum for T ~ e > but 
is nearly monotonic for e < 0. This distinguishes the 
model from transport through a single level. 

The inset of Fig. [I] shows the gate voltage dependence 
for different AB phases and Se = 0. The RG predicts A 
to decrease with increasing flux. For (j) — > tt the tunnel 
splitting is small and the level shift by A leads to a res- 
onance position of the linear conductance near e = —A. 
In contrast, for (f> — > 0, the tunnel splitting is large and 
?i w e which leads to a resonance position near e = 0. 
As a consequence we find that the position of the reso- 
nance is strongly influenced by the AB phase and reflects 
directly the tunnel splitting |A| together with the level 
renormalization A. 

Finally, we would like to comment on the case when 
the number of levels is given by N > 2. Generalizing 
the RG equation (||) to this case gives rise to an up- 
ward shift of all particle and hole excitations by approx- 
imately ~ AT, while only one level with an equally in- 
creased broadening remains approximately at the origi- 
nal position. This means that transport appears to be 
effectively controlled by single-level physics and may ex- 
plain recent experimentso'0 in the regime r ~ 6e where 
universal Kondo behavior of single-level dots has been 
observed. 

Summary. We have studied interaction and interfer- 
ence effects in quantum dots with two levels or two quan- 
tum dots with one level coupled to reservoirs. We found 
a new tunnel splitting that changes as a function of an 
applied magnetic flux and can be measured via the ab- 
sorption power. As function of the flux, the system can be 
tuned into an effective model showing Kondo physics. We 
expect important implications of our results for transport 
and spectroscopy experiments as well as for the theory 
of level statistics in quantum dots. 
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FIG. 1. Spectral density of the /i level for ei = T = 0, 
U = 10r, Se = 0.25, 0.5, 1, 2r (from left to right). Inset: 
5e = T, T = 0, 0.5, 1, 1.5r (from top to bottom). 



0.015 



0.01 

3 

< 
U 

0.005 



0, 



1 1 1 

<\ 

1 \ 0.2 


1 


1 1 


A 

\ 


— 8e/r=0.3 
-0.5 


- / v\ 

■ h ' s 




- 1.0 


1 . . 




-4 


4 


p\ i » \ 






/ \ / V \ 




— 5e/r=0.2 




-0.5 


/' „ A - — 
-4f A ^-^ ~- 


N 


- 1.0 


f'J V_ x ^ 

L ' ■ £ , 1 


1- 


■ — 

~" -T 







co/r 



FIG. 2. Absorption power (scaled, single peaks) vs. spec- 
tral density of the /i level ("shoulder" at the same position) 
for ei = -10r , U = 50r, and T = 0. In the inset, the 
spectral density is shown for ei = — 3.5r and U = 1QT. 




FIG. 3. Effect of a finite AB phase on the single parti- 
cle spectrum. The spectral density of the /i level is shown 
for ei = -1.6r, 5e = 0, U = 8.1T, and T = 0. Left in- 
set: Partial spectral density corresponding to the level C\. 
Same parameters as above, but with a finite level splitting 
Se = 0.08r. Right inset: /i spectral density for different 
broadening strengths Ti 7^ T2 of levels ci and C2. 
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FIG. 4. The linear conductance for 5e = 0, D — 50T, 
and U = oo. Main panel: T-dependence with = (dashed 
lines indicate the positive energies). Inset: ^-dependence for 

r = r. 



